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Abstract. The characteristic size of penumbral structures are still below the 
current resolution limit of modern solar telescopes. Though we have seen a sig- 
nificant progress in theoretical work over the last decades no tight constraints 
can be placed on the size of penumbral structures in order to favor models with 
relatively large and thick magnetic flux elements, just at or below the current 
resolution limit, or on the other hand, clusters of optically thin micro-structures. 
Based on a macroscopic 2-component inversion and the approach of polarized ra- 
diative transfer in stochastic media, we have estimated the characteristic length 
scale of the magnetic fluctuation in a sunspot penumbra from observed Stokes 
spectra. The results yield a coherent picture for the entire magnetic neutral 
line of the penumbra and indicate that the magnetic fluctuations have a typical 
length scale between 30 km and 70 km. 



Introduction 



Solar sunspots are the most prominent manifestation of concentrated magnetic 
fields on the solar surface and have been the subject of an exhaustive amount 
of research. Despite the fact that the photospheric structure and morphology 
is well studied and describ ed, and despite the advances in theoretical work (see 
e.g.[Th omas fc Weissl[2003 ) there is no generally accepted picture available which 
can explain the formation, structuring, and dynamics of sunspots. It is generally 
accepted that one of the keys to increase our overall understanding of sunspots 
lies in a better understanding of the many observable fine-scale features. The 
penumbra with its convectively driven small-scale filamentary structure provides 
an excellent test-bed for any model of magneto-convection. A wealth of observa- 
tions and modelling techniques gave us a good insight into the overall pattern and 
geometry of the magnetic field and plasma flows in the penumbra. The improve- 
ments in resolution and polarimetric sensitivity in the recent years has revealed 
an increasing number of small-scale variations within the penumbral magnetic 
field (strength and i nclination), velocity, and temperature structures (see e.g. 
Schlichenmaier 20021 ). Recent imaging observations have even revealed fibrils of 
150 to 180 km in width which contain subs tructures in the form of dark narrow 
cores, and which are less than 90 km in size (jScharmer et alj|2002l ). A small mag- 
netic structuring on sub-resolution scales (perpendicular and along the LOS) 
has also lon g been suggested and recognized by the ubiquitous Stokes profile 
asymmetry ( Solanki fc Montavon 1993 1 Schlichenmaier et al. 20021 ; iMiiller et al.1 
20021 : iTritschler et al.ll2004l : ISanchez Almeidall2005l ). Based on the specific sig- 
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natures of unresolved magnetic field structures in the Stokes profiles, we ex- 
ploit the diagnost i c cap a bilitie s of the meso-structured approach proposed by 
Carroll fc Staudd (|2003l . I2005T ): ICarroll fc Ko^l (120061 ). which is based on a 



stochastic formulation of the polarized radiative transfer. This new approach is 

applied to spectro-polarimetric observations of a sunspot penumbra which were 

subjec t to a macroscopic 2-component inversion bv lBellot Rubio. Balthasar. Sz ColladosI 
(2004|, hereafter RBC). Using this combination of macroscopic and mesoscopic 
inversion allows us to estimate the characteristic length scale of the underlying 
magnetic fluctuation. 

This paper is organized as follows: In Sect. 2 we give a brief overview of 
the basic concepts of the line formation in fluctuating and stochastic media with 
structures of finite length. In Sect. 3 we demonstrate that the statistical scatter- 
ing and absorption terms of the stochastic transfer equation directly determine 
the degree of the area and amplitude asymmetries of the resulting Stokes- U 
profiles, which can be used to estimate the typical length scale of the magnetic 
fluctuation. In Sect. 4 we describe the basic results of the 2-component inversion 
and the following application of the meso-structured approach. Finally, in Sect. 
5 we summarize our results. 



2. The Stochastic Transfer Equation for Polarized Light 

Any approach to radiative transfer modelling has to make assumptions about 
the underlying atmospheric structure. In atmospheres which are composed of 
small scale fibril structures - smaller than the resolution element - the natural 
question arises how these structures can be adequately accounted for in po- 
larized radiative transfer modelling? The two (magnetic) modelling approaches 
based on a macro-structured (conventio nal flux-tube model) or micro-structured 
(MISMA, ISanchez Almeida et al. 19961 ) description of the atmosphere cannot 



deal with a-priori unknown structures of finite spatial extent. To describe the 
polarized radiative transfer in an arbitrary fluctuating atmosphere, we pursue 
a rigorous probabilistic formulation of the underlying atmospheric structure. 
A detailed description of that stochastic approach and the derivation of the 
stocha s tic transport equation for polarized light, is given in ICarroll Staude 
(2005); [Carroll Kopj (|2006l ) . In this contribution we will only summarize the 



basic concepts of this approach which puts forw ard the idea of a more general 
meso-structured mag netic atmosphere (MESMA. ICarroll fc Kopll2006l ). 

A prerequisite for a stochastic modelling of structures with finite spatial ex- 
tent is to take into account the complete hierarchy of spatial correlations. If we 
assume that the spatial coherency or correlation of the individual atmospheric 
components along a given trajectory can be approximated by a Markov process, 
all the higher spatial correlations can be reduced to first-order correlation effects. 
We therefore introduce a random atmospheric vector B which comprises all rel- 
evant atmospheric parameters such as temperature, pressure, velocity, magnetic 
field strength, magnetic field inclination, etc., and specify the spatial correla- 
tion of the individual atmospheric components by the following discontinuous 
Markov process (the space-dependent, conditional-probability density function) 

p(B", s + As\ B\ s) = e -^ As S(B' - B") + [l - e -^ As ] p(B") , (1) 



The Magnetic Fluctuation in a Sunspot Penumbra 



3 



where j(B') is the fluctuation rate of the atmospheric structure B s and p(B") 
the spatial independent (stationary) probability density. The fluctuation rate 
can be expressed with the help of the correlation length I (characteristic length 
scale) to give 

7 (B s ) = r 1 (B s ). (2) 

This discontinuous Markov p rocess, also known und er the name Kubo- Anderson 
process, was already used by lFrisch Sz Frischl ( 19761 ) to describe inhomogeneous 



velocity fields in stellar and solar atmospheres. This process describes the spatial 
relationship in terms of the conditional probability density. p(B" , s + As \ B', s). 
It gives the probability of being in the state or atmospheric regime B" at the 
spatial position s + As, after having moved a short distance As from the known 
atmospheric state B' at s. The probability for staying in the initial regime B' , 
expressed by the Dirac delta function, decays exponentially with As while a 
sudden jump into another regime B s+ & s has an exponentially increasing prob- 
ability and is only weighted by the overall stationary probability density of the 
final state B" . This conditional probability density function thus describes the 
correlation between the two spatial positions s and s + As along a given trajec- 
tory (later the LOS). The degree of correlation is controlled by the fluctuation 
rate, 7, or the correlation length, I, respectively, which describe the mean path 
length or extent of the atmospheric structure. This generic stochastic model 
of the atmospheric structure allows us to formulate a differential equation for 
the depth-dependent atmospheric probability density function (PDF) p(B, s) 
from which we can derive a stochastic transp ort equation for polarized light (see 
ICarroll fc Stauddl2005l : ICarroll Ko P jl2006l . for a detailed derivation). For the 
particular application in this contribution, an appropriate model may be given 
by the discontinuous Kubo- Anderson process Eq. ([Tj) , which allows to derive a 
master-like transport equation for the mean conditional Stokes vector, which 
reads 



dY B 



K B Y B 



<h ~~ J1 ...„+jB+Jj(B)Y B ,p(B', S )dB'-J 1 (B')Y B p(B')dB', (3) 
where 7 is the modified fluctuation rate 

1(B) = t^tw, , (4) 
p(B,s)l(B) 

and Y B is the mean conditional Stokes vector which is defined as 

Y B (s) = J Ip(I,s\B,s) dl . (5) 

Please note, that p(B) and p(B,s) are in general not the same distributions, 
the latter is the spatial dependent pdf of the atmospheric state vector B which 
evolves from a given initial state along the ray path, whereas the station- 
ary pdf, p(B), is an a priori assigned attribute of t he atmosphere. The evo - 
lution of p(B,s) is g overned by a Master equation ICarroll &: Staudel |2005); 
Carroll &: KopJ (2006). If the atmosphere is structured on finite scales, p(B, s) 



eventually converges to the stationary distribution p(B), the equilibrium state. 
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Figure 1. Left: Area asymmetry vs correlation length. Right: Amplitude 
asymmetry vs correlation length. The area and amplitude asymmetries reveal 
a clear functional dependence on the correlation length. The dashed lines in 
both figures indicate the area and amplitude asymmetries which would origi- 
nate from a pure micro-structured atmosphere, calculated under the MISMA 
approximation. 



The stochastic transport equation for the mean conditional Stokes vector 
has a rather simple physical interpretation. Since Yb is conditioned on a particu- 
lar atmospheric regime B, and the transport equation (J3j) describes its statistical 
evolution through the atmosphere, four physical processes govern the transport 
of the mean conditional Stokes vector: the two processes of true absorption 
and thermal emission, and two processes that describe the statistical inflow and 
outflow of intensity to or from the regime B under consideration. In detail, the 
third term on the r.h.s. of Eq. ([3]) gives the amount of intensity or photons which 
may enter from all other atmospheric regimes into B and can therefore be con- 
sidered as an additional (statistical) source or scattering vector, while the fourth 
term describes the (statistical) loss or absorption of intensity due to transitions 
of photons from B to all other possible atmospherical regimes. The degree of 
statistical scattering and absorption is controlled by the correlation length I of 
the particular atmospheric regime. The observable of our problem — the expec- 
tation value of the Stokes vector at the top of the atmosphere — can easily be 
obtained from a final integration of the mean conditional Stokes vector over the 
entire atmospheric state space of B, 

(I(s))= jY B ( S )p(B,s)dB . (6) 



3. Stokes Profile Asymmetries 



Stokes- V profile asymmetri es can be attributed to the existence of gradients in 
veloc ity and magnetic fields (Landol fi Landi Degl'Innocent!ll996l ; Lopez Aristd 
2002). Thus the net circular polarization (NCP) and area asymmetry are 
good indicators for t he underlying magnetic inhomogeneity along the LOS. 
Carroll fc Kopj (120061 ) show that in a magnetic fluctuating or randomly or- 



ganized atmosphere the statistical scattering term in Eq. (j3J) is responsible for 
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Figure 2. Left: The averaged Stokes - V profile from the 2-component in- 
versio n (in the macroscopic limit) by iBellot Rubio, Balthasar. fe Colladosl 
(j2004l ). Right: The averaged Stokes-F profile after adjusting the correlation 
length to 40km (diamonds: observations; solid lines: fitted profiles). 



most of the NCP. To demonstrate this effect, we performed some simple model 
calculations for the Fe I line at 1564.8 nm where we have assumed a stochastic 2- 
component atmosphere, which consists of a field-free ensemble of structures with 
a net downflow of 1.5km/s, and a stationary magnetic ensemble of structures 
with a field strength of 700 G. To quantify the degree of asymmetry we used the 
usual definitions of the area asymmetry, 5 A, and the amplitude asymmetry, 8a, 

h V Wd\ max{y(A)}+min{y(A)} 
h\ V ( X )\ dX ' max{y(A)} - mm{V(X)} ' {> 

where s is the sign of the blue lobe of the Stokes- y profile. 

Figure Q] shows the Stokes- y profile area and amplitude asymmetry for 
Fe I 1564.8 nm as a function of the correlation length. Figure Q] nicely demon- 
strates, how the area and amplitude asymmetries depend on the correlation 
length (characteristic length scale) of the underlying atmospheric structures. 
The smaller the correlation lengths, or higher the fluctuation rates, the stronger 
the asymmetries. Please note, that only the correlation length has changed in 
that scenario. For very small correlation lengths the asymmetries asymptoti- 
cally reach the micro-turbulent limit (I — > 0). For a better comparison of the 
asymptotical behavior the true micro-turbulent limit was calculated under the 
MISMA approximation and is denoted by the dashed line in Fig. [TJ For very 
large correlation lengths and low fluctuation rates (the macrostructured limit) 
the asymmetries asymptotically approach zero. This is the expected result as 
the magnetic and non-magnetic structures in our model calculations possess no 
intrinsic gradients. This behavior is also a direct consequence of the stochas- 
tic transport equation (J3]) where the statistical scattering and absorption are 
inversely proportional to the correlation length I, and in the limit I — > oo the 
different atmospheric regimes finally decoupl e and no statistical scattering and 
absorption can occur. ICarroll fc~K opf (2006) show that the stochastic transfer 



equation presented here includes, in fact, the micro-structured (MISMA) ap- 
proximation as well as the macro-structured approach as limiting cases. Please 
note, that already for correlation lengths greater than 10 km significant devia- 
tions from the micro-turbulent approximation occur. As the inherent magnetic 
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fluctuation of the atm osphere has an impor tant effect on the degree of the NCP 
or area asymmetries ( Carroll fc Kopf 20061 ) . these asymmetries can be used in 



turn to estimate the underlying length scale of the magnetic fluctuation. 



4. Penumbral Structures from a Stochastic 2-Component Approach 

The starting point of our analysis are the results of a 2-component inversion 
of a sunspot penumbra from RBC. Their retrieved overall penumbral model is 
in agreement with a discontinuous mode l of the so-called uncombed penumbra 
proposed by ISolanki Sz Montavon ( 19931 ). with an almost horizontal flux-tube 



component embedded in a more vertical background field component. The in- 
version were performed by a special adaptation of the SIR code (see RBC for 
details) and returns the temperature stratification, the three components of the 
single-valued magnetic field vector, and the LOS velocity of the two components 
of the model, together with a single value for the macro-turbulent velocity. In 
addition, the stray light factor a and the filling factor of the flux-tube compo- 
nent are also determined. We do not go into the details of the analysis and 
results of the inversion of RBC, but we want to emphasize that the quality of 
the fitted profiles of the 2-component inversion seems to be able to grasp the es- 
sential character of the underlying penumbral atmosphere, despite the fact that 
the magnetic field and velocity structures of the two components are assumed 
to be neither depth dependent nor interlaced. The 2-component inversion of 
RBC represents therefore a typical macro-structured approach, because the two 
components are not interlaced along the LOS and, unlike the observed Stokes- F 
profiles, the resulting synthetic profiles therefore exhibit no NCP or area asym- 
metry. In order to allow for a more flexible meso-structured approach we have 
incorporated the model parameters of the 2-component atmospheres into our 
stochastic (meso-structured) transfer model. The underlying atmosphere in this 
stochastic scenario can now be considered as the composition of two different 
ensemble structures (the background and the flux-tube ensemble) but this time 
they are really interlaced. The mean length scale or extent of the individual 
ensemble structures is determined by the correlation length, which at the same 
time determines the rate of the fluctuation within the atmosphere and the abil- 
ity of the atmosphere to produce a NCP. For each individual average profile 
the correlation length of the underlying structures were adjusted until a best 
fit with the observation were reached. The profiles were averaged over a small 
portion of pixels (3 x 3). We put particular emphasis on profiles along the 
magnetic neutral line because of the preferred viewing angle onto the magnetic 
structures in this region. Figure [2] shows an averaged Stokes- F profile located in 
the magnetic neutral line. Despite the fact that the initial model parameters of 
the 2-component inversion are relatively good (Fig.[2l left), they cannot account 
for the area asymmetry of the observed profile. The observed Stokes- y profile 
of Fe I 1564.8 nm has a NCP of 4.9 mA and an area asymmetry of —17.6%. 
Since the infrared iron line is not very sensitive to velocity and magnetic-field 
gradients or fluctuations, due to its weakness and small range of formation, the 
observed profiles must be subject to a rather large fluctuation rate. After ad- 
justing the correlation length (with all other parameters kept fixed), the fit of 
the profiles could be considerably improved and the resulting NCP of 4.9 mA 
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and area asymmetry of —17.3% could be almost exactly reproduced. To obtain 
these large values of the Stokes- y asymmetries, with the given model parameters 
of the 2-component inversion, a correlation length as low as 40 km was needed. 
The analysis of a number of Stokes- y profiles along the magnetic neutral line 
gave similar results. For each of the analyzed Stokes- y profiles we could obtain a 
significantly improved fit compared to the macroscopic analysis of RBC, simply 
by varying the underlying correlation length. From our analysis we could infer 
that the typical length scale of the magnetic fluctuation in a penumbra must be 
between 30 and 70 km. 



5. Summary 



The st o chastic polarized radia tive transfer recently developed bv lCarroll Staude 



2005h: ICaTroll &; Kopfl (12 006) was applied to a meso-structured magnetic model 



(MESMA. ICarroll fe Kopf 20061 ) of a sunspot penumbra to estimate the char 



acteristic length scale of the magnetic fluctuation. The degree of fluctuation 
directly controls the amount of the NCP (or area asymmetry) of the resulting 
circular polarization signal. A maximum degree of asymmetry is always pro- 
duced by a micro-structured or micro-turbulent atmosphere, whereas the NCP 
vanishes for a macro-structured atmosphere. However, as it can be seen from 
Fig- [H these two cases only represent the very far limits of a more general 
meso-structuring that is valid for the broad range of finite length scales. Based 
on the results of 2-component analysis of a sunspot penumbra (RBC), we ap- 
plied a method which is able to account for the finite fluctuation of magnetic 
components along the LOS. With this approach, we were able to estimate the 
correlation length of the underlying magnetic structures, for a number of lo- 
cations along the magnetic neutral line. For the characteristic length scales of 
the magnetic fluctuation we found values between 30 and 70 km. These findings 
indicate that the magnetic field in the penumbra, although it may have on av- 
erage a typical uncombed-like structure, may possess a considerable amount of 
substructures. We conclude that the meso-structured analysis presented here, 
may provide a viable approach for studying small-scale and complex magnetic 
fields in the solar atmosphere. 
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